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Abstract. We present an estimate for the non-linear parameter tnl, which 
measures the non-gaussianity imprinted in the trispectrum of the comoving curvature 
perturbation, Q. Our estimate is valid throughout the inflationary era, until the slow- 
roll approximation breaks down, and takes into account the evolution of perturbations 
on superhorizon scales. We find that the non-gaussianity is always small if the 
field values at the end of inflation are negligible when compared to their values at 
horizon crossing. Under the same assumption, we show that in Nflation-type scenarios, 
where the potential is a sum of monomials, the non-gaussianity measured by tnl is 
independent of the couplings and initial conditions. 

Keywords: Inflation, Cosmological perturbation theory, Physics of the early 
universe. 
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1. Introduction 

Recent observations of the cosmic microwave background (CMB) have provided strong 
evidence [USE] that an era of inflation 01 El El El El HD] occurred in the early history of 
the universe. During inflation, the universe underwent a phase of accelerated expansion 
and the amplification of vacuum fluctuations in light fields generated a spectrum of 
density perturbations [HI [121 DEI EE] • This spectrum is conventionally parametrized in 
terms of the comoving curvature perturbation, (, on surfaces of uniform energy density. 

A central question in modern cosmology is to determine the physics which gave rise 
to inflation. It has recently become clear that a great deal of important information is 
encoded in the non-linearities of the curvature perturbation [151 EE El EH1 US]- The 
lower-order moments of ( are conventionally expressed in terms of the non-linearity 
parameters / NL $M EH E2] and r NL [SSI El ESI EH E7j , which respectively quantify the 
importance of the three- and four-point correlators of ( relative to the power spectrum. 
To date, primordial non-gaussian features in the CMB have yet to be detected and 
the current upper bounds on the non-linearity parameters are |/nl| i$ 10 2 [H EH EH] 
and |tnl| ^ 10 8 [27]. A major goal for the next generation of CMB experiments will 
be to improve on these bounds, and the Planck satellite is expected to be sensitive to 
|/nl| ~ 5 [21] and |tnl| ~ 560 [21]. Furthermore, it has recently been suggested that 
a future all-sky survey of the 21-cm background has the potential to probe down to 
|/nl| < 0.01 [30]. 

Given these forthcoming improvements in the observational constraints, there is 
a pressing need to obtain theoretical predictions for / NL and r NL in well-motivated 
inflationary models. In the single-field scenario driven by a canonically normalized 
scalar field, ( is conserved after horizon exit. This implies that non-gaussianities cannot 
be generated on scales larger than the horizon. Some non-gaussianity is imprinted by 
quantum mechanical interference at the epoch of horizon crossing [31] , but the level of 
this effect is small. In general, one expects /nl ~ t nl ^ rVlO for canonical single-field 
inflation [20j E2] , where r < 1 is the tensor-to-scalar ratidj. 

In multiple-field models, on the other hand, ( can evolve after horizon exit due to 
the presence of isocurvature modes [361 EH EH]- This allows additional non-gaussianity 
to be generated [261 EHJ HQ]. This possibility has been investigated by several authors, 

leading to estimates for / NL @H S21 E3 113 112 S3 SSI El SHI S3 EDI E3 EH E21 E3E11E3 

1561 157] 158] and tnl (23 EHJ E21 E3 [53] for a number of scenarios. In particular, Vernizzi 
& Wands [H] have determined the form of /nl i n a two-field inflationary model with 
a separable potential. They find that |/nl| is typically small in the so-called "horizon- 
crossing approximation" , where the fields roll to their minima after inflation with values 
that are negligible compared to those at horizon crossing. The same conclusion was 
reached by Alabidi & Lyth [27], and by Kim & Liddle [B] (see also Sasaki, Valiviita & 
Wands |53j), who used the horizon-crossing approximation to compute the non-linearity 

| Much larger effects are possible in models with non-canonical scalar fields [331 1311 [35] , but we do not 
consider such scenarios in the present work. 
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parameters for multi-field models with monomial potentials. Recently, Battefeld & 
Easther [57] have extended the analysis of Vernizzi & Wands to an arbitrary number of 
scalar fields and have found the corresponding form of the non-linearity parameter /nl 
without assuming the horizon- crossing approximation. 

In the present paper, we derive the first complete, tree-level estimate of the 
trispectrum non-linearity parameter, Tnl, for inflationary models driven by multiple, 
uncoupled scalar fields. In §2] we briefly describe the non-linear 5N formalism, and 
outline how it can be used to compute tnl- In particular, we present the complete, 
tree-level SN expression for t^l in §2.21 In $31 we discuss how the derivatives of N 
are computed in a scenario with an arbitrary number of uncoupled scalar fields and 
we assemble these components in §|4] to derive an exact analytic formula for the non- 
linearity parameter. Our result is valid in the slow-roll approximation and accounts for 
the evolution of perturbations on superhorizon scales. In §j5j we estimate the magnitude 
of tnl in scenarios driven by monomial potentials. We conclude with a discussion in §|6j 

Throughout this paper we work in natural units where the Planck mass, defined by 
Mp 2 = 8nG, is set to unity. 

2. The non-linearity parameter of the trispectrum 

2.1. The 5N formalism 

The non-linear 5N formalism provides a powerful method for calculating non- 
gaussianities in <Jjj. It is valid under the assumption that causally disconnected regions 
of the universe evolve locally like an unperturbed universe in which the pressure and 
density are allowed to vary [62j EH El] • This leads to the relation ( = 5N, where N((f>, p) 
represents the number of e-folds of expansion between some flat primordial slice, *, at 
time on which the fields {4> a } have assigned values, and a final uniform-density surface, 
c, at time t c on which the energy density has the specified value p [59j [60l [36l [39j [26]. 

If slow-roll inflation is valid on the initial flat slice, the curvature perturbation can 
be written as a power series in the field perturbations evaluated at that time |39j : 




where = dN/d(p aif1 iV a/ g* = d 2 N '/ 'dipa^dcj)^, and a subscript * denotes evaluation 
at £*. Note that although is a derivative with respect to the fields on *, it depends 
on the values of the fields on both * and c. This is also true for the higher derivatives 



§ The linear formulation was first developed by Starobinsky [59] and subsequently employed by Sasaki 
& Stewart [60] to calculate the perturbation spectrum for multiple-field inflation. It was recently 
generalized to the non- linear regime by Lyth, Malik & Sasaki [36]. (See also Langlois & Vernizzi [61] , 
and Rigopoulos & Shellard [IB] for alternative, but equivalent, formalisms). 





(1) 



in Eq. (□). 
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We consider inflation where M scalar fields contribute to the energy density of the 
universe through a separable potential]]] 

W(<P) = Y J VMe)- (2) 

a 

The scalar field dynamics can be determined in terms of the slow-roll parameters 

c i w ta w, aaa 

The slow-roll approximation is valid when the e a and \rj a \ <C 1. Another important 
quantity is e = —H/H 2 , which measures how far the inflationary era lies from de Sitter 
space, with e = for exact de Sitter and e < 1 for inflation. When the slow-roll 
approximation applies, one can express e in terms of the individual e a such that 

e = J>. (6) 



Q 



In this regime the equations of motion for the scalar fields reduce to 

3H<p a = -W >a . (7) 

Since these are first-order equations, the phase space of the system is M- dimensional. 

Finally, the number of e- folds of expansion N, whose derivatives appear in Eq. (pQ), 
can be written ClS db quadrature: 

-fx: < s > 

2.2. The four-point correlator 

The power spectrum of ( is formed by combining two copies of Eq. ([T]) and taking the 
expectation value in the quantum vacuum state. If only the lowest derivatives of iV are 
retained, it follows that 

(C(ki)C(ka)) = (27r) 3 5(k 1 + k 2 )P,(A; 1 )^iV i L, (9) 

a 

where P*(ki) is the power spectrum of a massless scalar field evaluated on * and is 
related to the "dimensionless" power spectrum "P* by 

k 3 P,(k) Hi 

'* [ ) 27T 2 47T 2 ' 1 ' 

The four-point correlator is derived in a similar fashion by forming a product of four 
copies of Eq. ([1]) and taking expectation values. In general, the trispectrum is defined by 
(C(ki)C(k 2 )C(k 3 )C(k 4 )) = (2Tc) 3 5(Y.i ki)T(ki, k 2 , k 3 , k 4 ), where only the connected part 

|| In this paper we will not use the summation convention on scalar field indices {a, (3, 7, . . .}. 
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of the correlator is considered [2"3"l [21]. Observational constraints on the trispectrum 
are conventionally given in terms of a non-linearity parameter, t^l, which is defined by 
[251 ESI EU 

T(k 1 ,k 2 ,k 3 ,k 4 ) = ^T NL ( y J2 N l*)\ P *( k i) P *( k 2)P*(h4) + 23 permutations], (11) 

a 

where k^ = kj + k,. In general, tnl is a function of the four momenta {ki, k 2 , k 3 , k 4 } 
which parametrize the trispectrum. It is convenient to rewrite Eq. (fTTj) in such a way 
that the momentum-dependence is made explicit. If we define a function T such that 

we find that [32] 

^=^N L (E^) 3 ^ T - (13) 

After taking the product of four copies of Eq. (CQ) and evaluating the resulting 
convolutions, it can be shown that there are four distinct terms which contribute to 
and in which unconstrained integrals over momenta do not appear [32] . We will refer to 
such terms as the tree-level terms. The leading tree-level term was computed explicitly 
in Ref. [32] and found to have the form 

^4 



At, 



(i) 

NL 



r(E*N^y (14) 

where T\ = ^ perms -M 4 (ki, k 2 , k 3 , k 4 ), and M.^ is a momentum-dependent form-factor 
given by [32J 



, kjkj W 24 
'P P K 

I hi 



Z^ 2 • Z 34 3 

— r2 h 2k 2 • Z 34 + -<Ji 2 (j 34 



2 k 



2 °34 
31 



k x ■ k 2rTr t k\k\ 



W 12 A + 



2 + 6 



Av 



The three momentum-dependent scalars a^, Wij and Wg mn are defined by 



a 



ij 



kj ■ kj + kj . 



W i:j = 1 + 



w £mn = 1 



/c^ ~\~ kj 2 /c^ /c j 



kt 



kl ' 



ki + fc m + k n 2(kgk m + /c£& n + fc m A; n ) Qkgk m k n 



Af 



A:? 



(15) 

(16) 
(17) 

(18) 



and Zij satisfies 



<7jjkj — o-jikj. The total scalar momentum is written k t = 
k\ + k 2 + k 3 + A; 4 . 

It was further shown in Ref. [32] that |Atj^l| < r/50, where r < 1 is the tensor-to- 
scalar ratio. This term is therefore small. 
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The remaining three contributions were also presented in Ref. [32J for single-field 
inflation. For a multiple-field scenario, they are given by 

NL= (E^.) 2 ' ' 

A (3) Ea^ 7 ^,«/3*^,a7*^,/3*^,7* , . 

At-stt = — 3 , (20) 

and 

A (4) = 2 E Q/ 3 7 N^^N^N^N^ ^ kf (21) 
(EA) 3 T > 

respectively, where K, is defined by 

K = ^ H ^M 3 (k 12 ,h,h) (22) 



T kin 

perms x ^ 



and 



k 2 k 2 1 PP 

„2 /( A 2"'3 , / 3 i ^2^3, 



M 3 (h, h, k 2 ) = -A;^ - 4^ + + -^(fe - fc 3 ). (23) 

t t 

Eq. ( J20l was derived previously by Alabidi & Lyth [27] and Lyth [51]. An expression 
equivalent to Eq. (|2ip was given by Sasaki, Valiviita & Wands [53] in the curvaton 
scenario. 

In general, there are an infinite number of other contributions to tnl which involve 
'loops' in the sense of unconstrained integrals over momenta [121 [51]. Such loops 
are suppressed relative to the tree-level by powers of the spectrum, V* < 1CT 10 . In 
principle, however, these loop integrals may be sufficiently large for such terms to become 
important and it is not known whether they are always negligible in comparison with 
the tree-level contributions. On the other hand, it was recently shown by Sloth that the 
contribution from a specific one-loop correction to the power spectrum is only about 1% 
that of the tree-level effect [65J. In view of this, we will assume that the contribution of 
these corrections is indeed negligible. 



3. The derivative of N 



Our aim is to derive expressions for the contributions (fl4"]l — (I2TT) in terms of the 
parameters ([3]) ([5]) for a general multi-field, slow-roll inflationary scenario when the 
potential is separable. The central problem is to compute the derivatives of iV which 
appear in Eq. (CEJ). In the horizon- crossing approximation, one ignores any contribution 
from c, in which case the derivatives can be written entirely in terms of quantities 
evaluated on * [66l SH (27J [44] . The method of Vernizzi & Wands [43] evades this 
assumption and was recently extended to an arbitrary number of uncoupled fields by 
Battefeld & Easther [57]. An alternative approach was followed by Lyth & Rodriguez 
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and by Alabidi [45] . but in this latter case it is necessary to explicitly integrate 
Eq. (J7j) between * and c. 

To proceed, let us consider a point in phase space that is specified by the values of 
the fields {<p a }- This point lies on a particular trajectory However, changing the values 
of the fields will move the point and consequently will also alter the trajectory. We 
need to compute the variation of iV under this change of trajectories. It follows from 
differentiation of Eq. (jSJ) that if the field values on * are changed by {d(fi a *}, the total 
differential of N is given by 



diV 



Va 



W/3 



<90/3c 



d0„*, 



(24) 



where subscripts * and c denote quantities evaluated on the initial flat hypersurface and 
the final comoving hypersurface, respectively. 

The values <ft ac are completely determined by the trajectory on which they lie. 
Since we are assuming that slow-roll inflation is still occurring at the final comoving 
hypersurface, the energy density p on c is dominated by W c , which is constant over the 
entire surface. Using this property, Battefeld and Easther show that the derivatives of 
(f) ac with respect to 0/3* are given by [5?j 

W c H~_ ( ep c 

a /3 




Wm y epm , . ( 25 ) 

In obtaining the derivatives of N, it will prove convenient to introduce another set 
of parameters, 

E a p = — — 5 aj 3, (26) 



Y, 



-7c 



1 - 



Vic 



(27) 



and 



X a 



(28) 



oc 

which are defined in terms of the values of the slow-roll parameters on c. We will also 
need expressions for the variations of the slow-roll parameters on c with respect to a 
field on *. In particular, we will make use of the following relationships: 



d(f>x* 
de c 

<90A* 



)E 



Xa j 



(29) 
(30) 



% This formula was derived in Ref. [57] given a specific choice of conserved quantities which label the 
allowed inflationary trajectories. In general, one can always rearrange the conserved quantities by 
taking arbitrary linear combinations and we have verified that Eq. (|25| applies for this more general 
choice. 
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and 

dVac V2 W c 





y/2 W c 








V2 W c 







, T/ £acE\ a . (31) 

These follow after differentiation of Eqs. (J3j) (J3D and substitution of Eqs. fT25]) - (l28]) . 
Finally, taking the derivatives of Eqs. (1261) and (|2"T|) with respect to <fi^ yields 

X^E^ + (Y 1 - e 7C ) E ^Ys\ , (33) 

where we have employed Eqs. fl28]) - fl3Tj) . 

The first two derivatives of iV were calculated in Refs. [I3j EZ] in terms of a 
parameter, Z ac , defined by 

Z ac = ^W c - V ac . (34) 

This quantity is evaluated purely on the final comoving hypersurface c and measures 
the importance of any isocurvature effects in the generation of non-gaussian statistics. 
It vanishes identically in single-field inflation and, more generally, ( is conserved on 
super-horizon scales whenever Z ac = 0. In this case, the derivatives of N are controlled 
entirely by the potentials and the slow- roll parameters evaluated on *. On the other 
hand, when more than one scalar field contributes to the energy density during inflation, 
( can evolve on large scales due to the contribution from Z ac . 

By employing Eqs. (124"1) and (1251) . it can be shown that the derivative of N with 
respect to a field value on * is given by 

N a * = (35) 

where 

V a:¥ -|- Z ac /nr>\ 

u a EE — — . (36) 

The second derivative of N is then derived by differentiating Eq. ( 1351) : 

It follows that N only has a non-trivial, mixed second derivative if there is a non- 
negligible isocurvature contribution. 

To compute the non-linearity parameter tnl, we require the Taylor expansion of iV 
up to and including the third derivative. We find after some algebra and use of Eqs. 
and (133|) that differentiation of Eq. (1371) yields 

^,a/3-y* = Oa/3-y , U a — a p- 77T~^~ 0<*7 



1 d^Z ac . . 
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where 5 a ^ is unity when a = (3 = 7 and is zero otherwise. 

Finally, we must relate the derivatives of Z ac to the slow-roll parameters. After 
differentiating Eq. ( 1341) twice with respect to the field values on * and making use of 
Eqs. (|26]) and fl2SH2T), we find that 

9Z - ^ W -A^ (39) 



* v (3* 

d 2 Z ac rjp* dZ ac 2W* 



where we have defined the totally symmetric matrices ^4 Q) g [331 EZ] and *4 Q/ 3 7 such that 



= — 77^7 2^ Y^asEpS, (41) 



<5 

and 



* s ^ 

+ 2J — e 5c) [Ea^E^sEys + Ef3 K E y sE a s + E 1K E a $Ep$] >. 



(42) 



Since iV^* and N )0l p 7 * are derived from the partial differentiation of a scalar 
quantity, they are symmetric under the exchanges a ±=> [3 and a ±=> [3 ±=> 7, respectively. 
Indeed, this symmetry is manifest since A a p and A a p-y are explicitly symmetric under 
an interchange of all indices. 

4. An analytic formula for Tnl 

Eqs. (I35D — (I38|) represent expressions for the derivatives of N given entirely in terms of 
the values of the potentials and slow- roll parameters ©-([5]) evaluated on c and *. It 
is now possible to derive an analytic expression for the non-linearity parameter tnl in 
terms of these parameters by evaluating the corresponding expressions for the quantities 
A Tnl, At N l and At N l defined in Eqs. (fT9|) - (!2T|) . We begin by writing 

Atnl= (*>) 

(£.£) ' 

where r e {1,2,3,4}. Eq. (fT4"j) implies that t*- 1 ** is given by 

and, by employing Eqs. (|35|) . (!37j) - (l38|) . (j4"Tl) and (f42|) . we find that the remaining 
are given by 



(45) 
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■aP 



a x ' a.p 

+ J] p ^ A^Ary, (46) 

a/?7 



^a*^/3*^7* 



and 

r 



8£A 3 



6 0-2 I, V'a* ?a*J '' a * J Z^O f 2 , 'la*^aP 

a zfc a* \ fc a* / a ^ z fc a* fc /3* 



£0*6/3*^7* 

respectively. 

Eqs. (l43j) - (H71) are valid at any time during slow-roll inflation and represent the 
complete, tree-level non-gaussianity imprinted in the trispectrum for an inflationary 
scenario driven by multiple, uncoupled fields. 



5. The horizon-crossing approximation 

In this Section, we compute the non-gaussianity in tnl for the class of inflationary 
models where the potential W(4>) is given by a sum of monomials: 

W(0) = X^"- (48) 

a 

This class includes the Nflation [63 EH [69] and assisted inflationary [Hfl EH E21 E31 E3! 
scenarios. Nflation [67J is a realization of assisted quadratic inflation within the context 
of string theory compactifications. The scalar fields that drive inflation arise as axionic 
degrees of freedom associated with shift symmetries in the compact manifold. Although 
the axion potentials are sinusoidal, they can be approximated by a quadratic form |A a 0^ 
when the fields are in the vicinity of their minima. 

In many inflationary models, such as Nflation and assisted inflation, the scalar 
fields settle into their minima during the final stages of accelerated expansion. In this 
case, the end-point (on the final surface c) of the integral jS]) is unimportant. This 
implies that the number of e- folds, N, becomes independent of the field values cf) ac . 
Moreover, since the parameter Z ac quantifies the effect of including this end-point on 
the evolution of the curvature perturbation, it too should be negligible. Ignoring such 
contributions is known as the "horizon- crossing" approximation and is equivalent to 
specifying Z ac m A a p ~ A a ^, ~ in Eqs. (H5l) - (|47|) . This implies that all quantities in 
the final expression for ( are evaluated at the time of horizon exit on the surface *. 

The horizon-crossing approximation has been employed previously to argue that 
the non-gaussianity imprinted in the bispectrum is always small in multi-field inflation 
[271 PfOl HI] . Alabidi & Lyth [27] have also used it to estimate the contribution t^ to the 
non-linearity parameter tnl- However, they did not employ the full tree-level formulae 
given by Eqs. (T44j)-(j47j). We will include all these contributions in what follows and 
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assume that the horizon- crossing approximation is valid. We will further assume that 
the coupling constants, A a , and powers, T a , in the potential (|48l) are arbitrary, subject 
only to the condition that a slow-roll phase of inflation can be realised. It then follows 
that the number of e-folds from the surface * to c is given by 



N 



( 49 ) 



We begin by considering the contribution from t^\ Since we are assuming that 
A a /3 ~ Aap-y ~ 0, it follows from Eq. (H9|) that t^ = 0. Hence, this term does not 
contribute to the non-linearity parameter tnl- On the other hand, the remaining terms 
dMD OS]) yield 

r T± r 2 ( v-^ A a 




where 



is the tensor-to-scalar ratio and Vq = V* J2 a ^ 2 a * * s the power spectrum of (. Eq. 
( 150]) depends explicitly on the couplings, A a , and the initial conditions, which are 
expressed through the field values on the initial slice *. Hence, the non-gaussianity 
measured by the trispectrum is model-dependent. 

Further conclusions may be drawn if the powers T a in the potential (148p take the 
same value, T a = T, as in Nflation for example. In this case, the tensor-to-scalar ratio is 
directly related to the number of e-folds such that r = 4T/N [27] . As a result, Eq. (1301) 
may be further simplified and we find that 

+ * < 52) 

Note that any dependence on the coupling parameters and the initial conditions of the 
fields has now disappeared. 

The strongest upper bound on r is imposed by the WMAP data [1], which implies 
that r < 0.55 at 95% confidence. The contribution for t^ is known to give a contribution 
to |Atnl| that is bounded from above by 0.06r [32]. Even assuming that r saturates the 
WMAP upper bound, this gives no more than |At^l| ~ 1/30, which is too small to be 
observable. The contribution from the other terms in Eq. ( l52l) will be maximized when /C 
is maximized and this is expected to occur on equilateral configurations. In this case we 
find that /C ~ —1.7. Hence, the WMAP bound on r implies that in the equilateral limit 
the second term in Eq. (|52|) can not exceed (3 x 10~ 3 )/r and this is also unobservably 
small. Although it is likely that some memory of the initial conditions will remain in the 
isocurvature-driven contributions A a p and A a ^, the considerations of Ref. [57] suggest 
that these too will be small, with the most significant contribution arising when the 
fields have a narrow spectrum of masses (or couplings). We conclude, therefore, that 
the non-gaussianity from assisted inflationary scenarios based on monomial potentials 
with equal powers is unlikely to be visible in the trispectrum. 
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6. Conclusions 

In this paper, we have employed the 5N formalism to derive the first complete estimate 
at tree-level for the trispectrum non-linearity parameter tnl in uncoupled, multi-field 
inflation. Our result assumes that the slow-roll approximation is valid between horizon 
exit, and the time of interest, which we have taken to be before the end of inflation. 
We have considered an arbitrary number of fields and accounted for the influence of 
isocurvature perturbations in generating non-gaussianities on superhorizon scales. 

We have employed our expression within the context of the horizon-crossing 
approximation to compute the non-gaussianity imprinted in the trispectrum by an epoch 
of inflation driven by a sum of arbitrary, monomial potentials with arbitrary coupling 
constants. In general, it was found that the non-linearity parameter tnl is model- 
dependent, since it is sensitive to the initial conditions and the values of the coupling 
parameters. However, when the powers in the monomials take the same value, the non- 
linearity parameter becomes independent of the initial conditions and couplings. In this 
class of models, we are able to estimate an upper bound on r NL which is determined by 
the tensor-to-scalar ratio. We find that the non-gaussianity is always small and will not 
be observable to upcoming CMB experiments, such as the Planck satellite. Alternative 
scenarios exist in which the non-gaussianity may be larger, such as the curvaton, for 
which a comparable prediction for the trispectrum exists [53]. A future measurement 
of the parameter tnl offers the possibility of a clean experimental comparison between 
these scenarios. 

Note added: After this paper was submitted, a similar work [75] appeared in which 
the authors use the 5N formalism to compute an expression for t^l comparable to 
Eqs. (|T9l)-(l2TT). and give an explicit slow-roll expression comparable to Eqs. (H6l)-( )47j) in 
the single-field case, when the bispectrum is zero and the horizon-crossing approximation 
applies. We have verified that in the appropriate special cases the results of Ref. [75] 
exactly agree with those presented in this paper. 
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